Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN style file vl.4) 



Cluster versus POTENT Density and Velocity Fields: 
Cluster Biasing and Omega 

E. Branchinii'2, I. Zehavi3'^ M. Plionis^'^ & A. Dekel^ 

^Department of Physics, University of Durham, South Road, Durham DHl 3LE, UK 

^ Kapteyn Institute, University of Groningen, Landleven 12, 9700 AY, Groningen, the Netherlands 

^ NASA/Fermilab Astrophysics Group, Fermi National Accelerator Laboratory, Box 500, Batavia, IL 60510-0500, U.S. 

"^Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

^National Observatory of Athens, Lofos Nimfon, Thesio, 18110 Athens, Greece 

^ SISSA - International School for Advanced Studies, via Beirut 2-4, 1-34013 Trieste, Italy 

1 February 2008 



ABSTRACT 

The density and velocity fields as extracted from the AbcU/ACO clusters are 
compared to the corresponding fields recovered by the POTENT method from 
the Mark III peculiar velocities of galaxies. In order to minimize non-linear 
effects and to deal with ill-sampled regions we smooth both fields using a 
Gaussian window with radii ranging between 12 - 20h-^Mpc. The density 
and velocity fields within 70/i"^Mpc exhibit similarities, qualitatively consis- 
tent with gravitational instability theory and a linear biasing relation between 
clusters and mass. The random and systematic errors are evaluated with the 
help of mock catalogs. Quantitative comparisons within a volume containing 
~ 12 independent samples yield f3c = il'^'^ /be — 0.22 ± 0.08, where be is the 
cluster biasing parameter at 15 h~^Mpc. If be ^ 4.5, as indicated by the cluster 
correlation function, our result is consistent with f7 ^ 1. 

Key words: Cosmology: theory - galaxies: clustering, large-scale structure, 
large-scale dynamics. 



1 INTRODUCTION 

The basic hypothesis underlying the study of large-scale structure is that it grew out of initial fluctuations via 
gravitational instability (GI). In the linear regime, this theory predicts a relation between the peculiar velocity and 
density fluctuation fields, V-U = —f{Q)S, with — OP'^' . From observations we can deduce the density field of 

galaxies or clusters rather than the density field of the underlying matter distribution. One then needs to assume 
a relation between the galaxy or cluster fluctuation field and that of the mass. A first order approximation is that 
of a linear "biasing" relation (hereafter LB) in which the two fields, smoothed on the same scale, obey the relation 
So — boS. Thus, GI-I-LB boil down to a simple relation between observables, 

V v = PoSo , i3o = n°-Vbo . (1) 

The density field of the extragalactic objects can be derived from a whole-sky redshift survey, while the velocity 
divergence can be reconstructed from a sample of redshifts and distances inferred by TuUy-Fisher-like distance in- 
dicators. Therefore, combining these data allow a measure of f3o, which, subject to some a priori knowledge of the 
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biasing parameter bo, provide constraints on the cosmological density parameter 57. A related analysis, invoking the 
integral of equation (|l|), can be performed using velocities rather than densities. 

The efforts to measure (3 from various data sets using different methods are reviewed in e.g., Dekel (1994, 
1997), Strauss & Willick (1995). The most reliable density-density analysis, incorporating certain mildly-nonlinear 
corrections, is the recent comparison of the IRAS 1.2 Jy redshift survey and the Mark III catalog of peculiar velocities 
yielding, at Gaussian smoothing of 12 /i~^Mpc, f3iRAS = 0.89 ± 0.12 (Sigad et al. 1998, PI98; replacing an analysis of 
earlier data by Dekel et al. 1993). An analysis of optical galaxies has provided a somewhat lower value for jSopt (Hudson 
et al. 1995), in accordance with the expected higher biasing parameter for early-type galaxies as demonstrated by 
their stronger clustering (cf. Lahav, Nemiroff & Piran 1990). Recent velocity- velocity comparisons typically yield 
values of Piras — 0.5 — 0.6 ± 0.1 (Willick et al. 1997b and references therein, Willick & Strauss 1998; Davis, Nusser 
& Willick 1996, da Costa et al. 1998, Branchini et al. 1999). The main source of uncertainty in the interpretation of 
the (3 estimates arises from our ignorance concerning the biasing relations. Fortunately, we do have a handle on the 
relative biasing parameters, based, for example, on the relative amplitudes of the correlation functions of the different 
types of objects, which should scale like fe^. Since different classes of extragalactic objects are assumed to trace the 
same velocity field, one can hope to tighten the constraints on Q. by deriving /3 for several different types of objects. 

Clusters of galaxies are promising candidates for this purpose because they are well-defined objects and are 
sampled quite uniformly to large distances, much larger than the available galaxy peculiar velocity samples. The use 
of the cluster distribution to probe the large-scale dynamics has been mainly restricted to dipole analyses, where the 
predicted velocity at the Local Group (LG) is compared to its observed motion relative to the CMB frame (Scaramella, 
Vettolani & Zamorani 1991; Plionis & Valdarnini 1991, PV91; Plionis & Kolokotronis 1998). It has been found that 
the directions of the two dipoles converge when using a large enough sample of clusters (> 150 /i~^Mpc), as expected 
from the assumed global homogeneity of the cosmological model (and contrary to the finding of Lauer and Postman 
1994, based on their attempt to directly measure peculiar velocities for clusters). Once the cluster distribution is 
properly corrected from redshift to real space, the corresponding value of (3 derived from the dipole is /3c — 0.21 
(Branchini & Plionis 1995, 1996, BP96; Scaramella 1995a; Branchini, Plionis & Sciama 1996). This estimate is higher 
than the value derived without this correction (Scaramella et al. 1991; PV91), and is consistent with ~ 1 for cluster 
biasing parameter of 6c ~ 4 — 5 as indicated by the cluster correlation analyses. The validity of the LB assumption for 
clusters might be questioned. In particular, such a large biasing parameter cannot follow the linear relation in deep 
underdensities. However, because of their low number density, clusters trace the underlying mass density field with 
a large inherent smoothing scale set by their mean separation. This has the effect of decreasing the density contrast 
and restoring the plausibility of the LB hypothesis over a large fraction of the volume sampled. 

It turns out that the bulk motion, as predicted from the cluster distribution with (3c ~ 0.2 inside a sphere of 
radius ~ 50/!.~^Mpc about the LG, is consistent with that derived directly from galaxy peculiar velocities (see Dekel 
1997, Dekel et al. 1998; Giovanelli et al. 1996, 1998a, 1998b). However, the estimate of /3c from the dipole at one 
point, or from the bulk flow, naturally suffers from severe cosmic scatter {e.g., Juskiewicz, Vittorio & Wyse 1990). The 
cosmic scatter can be reduced if the comparison is made at several independent points. Branchini (1995) and Plionis 
(1995) have attempted to compare predicted velocities from the cluster distribution to observed peculiar velocities of 
groups and clusters from Tormen et al. (1993), Hudson (1994), and Giovanelli et al. (1997), obtaining again /3c ~ 0.2. 
These analyses, however, are of limited validity since they compare smoothed and unsmoothed velocities. 

The purpose of this work is to measure /3c by comparing the Abell/ACO cluster distribution and the galaxy 
peculiar velocities of the comprehensive Mark HI catalog as analyzed by POTENT. The comparison is done alterna- 
tively at the density-density level and at the velocity-velocity level, and involves a careful error analysis. In § |^ we 
summarize the Mark HI data, the POTENT method and the associated errors. In § ^ we describe the reconstruction 
of the cluster density and velocity fields and the various sources of error. In § ^ we perform a quantitative comparisons 
of the cluster and POTENT fields, in order to determine /3c. We conclude our results in § |^. 
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2 POTENT RECONSTRUCTION FROM PECULIAR VELOCITIES 

The POTENT procedure recovers the underlying mass-density fluctuation field from a whole-sky sample of observed 
radial peculiar velocities. The steps involved are: 

(a) preparing the data for POTENT analysis, including grouping and correcting for Malmquist bias, 

(b) smoothing the peculiar velocities into a uniformly-smoothed radial velocity field with minimum bias, 

(c) applying the ansatz of gravitating potential flow to recover the potential and three-dimensional velocity field, 

and 

(d) deriving the underlying density field by an approximation to Gl in the mildly-nonlinear regime. 

The POTENT method, which grew out of the original method of Dekel, Bertschinger & Faber (1990, DBF), is 
described in detail in Dekel et al. (1998, D98) and is reviewed in the context of other methods by Dekel (1997, 1998). 
Farther improvements since DBF have been introduced which we use in the present analysis. They are discussed in 
detail by Sigad et al. (1998) 

We use the Mark 111 catalog of peculiar velocities (Willick et al. 1995, 1996, 1997a), which is a careful compilation 
of several data sets consisting of ~ 3000 spiral and elliptical galaxies. The non-trivial procedure of merging the data 
sets accounts for differences in the selection criteria, the quantities measured, the method of measurement and the TF 
calibration techniques. The data per galaxy consist of a redshift z and a "forward" TF (or Dn—a) inferred distance, d. 
The radial peculiar velocity is then u = cz — d. This sample enables a reasonable recovery of the smoothed dynamical 
fields in a sphere of radius ~50ft.~^Mpc about the Local Group, extending to ~ 707i~^Mpc in some well-sampled 
regions. 

The POTENT method is evaluated using mock catalogs. The mock catalogs and the underlying A''-body sim- 
ulation are described in detail in Kolatt et al. (1996, K96). Here we only stress that a special effort was made to 
generate a simulation that mimic the actual large-scale structure in the real universe, in order to take into account 
any possible dependence of the errors on the signal. 

2.1 Errors in the POTENT Reconstruction 

D98 and PI98 demonstrate how well POTENT can do with ideal data of dense and uniform sampling and no distance 
errors. The reconstructed density field, from input that consisted of the exact, G12-smoothed radial velocities, is 
compared with the true G12 density field of the simulation. The comparison is done at grid points of spacing 
5 h~^Mpc inside a volume of effective radius 40 /i~^Mpc. No bias is introduced by the POTENT procedure itself and 
they find a small scatter of 2.5% that refiects the accumulating effects of small deviations fi-om potential flow, scatter 
in the non-linear approximation and numerical errors. 

Using the mock catalogs described by K96, we want to check and quantify how well the POTENT reconstruction 
method works on our sparse and noisy data. Our goal is to eventually compare the POTENT fields to the density 
and velocity fields obtained from the distribution of clusters. Since the clusters are sparse tracers of the mass, we 
need to explore also smoothing radii larger than the G12 (commonly used in POTENT applications), and we check 
the G15 and G20 cases as well. The errors due to sparse sampling and nonlinear effects are expected to be smaller 
for the larger smoothing scales, while the sampling-gradient bias may increase. 

For each smoothing radius, we execute the POTENT algorithm on each of the 20 noisy mock realizations of the 
Mark 111 catalog, recovering 20 corresponding density and velocity fields. We will later consider the individual fields 
as well as the mean fields averaged over the mock catalogs. The error in the POTENT density field at each point 
in space, as^, is taken to be the rms difference over the realizations between Sp and the true density field of the 
simulation smoothed on the same scale. The errors on the smoothed velocity fields, ay^ , are estimated by a similax 
procedure from the Y supergalactic component of the mock velocity field. We evaluate the density and the velocity 
fields and their errors at the points of a Cartesian grid with 5 7i~'^Mpc spacing. In the well-sampled regions, with the 
G15 smoothing, the errors for the density are typically as^ «0. 1 — 0.3 and cr^p « 50 — 250kms~^, but they are much 
larger in certain regions at large distances. 
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The error estimates a^^ and a^^ are two of the criteria used to exclude the noisy regions from the comparison 
with the clusters. The third one is the distance from the 4-th neighbouring object in the Mark III catalog, i?4, 
which provide us with a measure of the poor sampling in the parent velocity catalog. Two more cuts have been 
applied on the cluster density and velocity fields, using the errors tj^^ and os^ obtained from the mock catalogs 
analysis in § 3.3.2 . Furthermore, we only consider objects within R = 70 h^^Mpc and outside the Zone of Avoidance 

> 20°). Our last constraint is on the misalignment angle between the cluster and the POTENT velocity vectors, 
A6', which we impose to be smaller than 45°. The reason for such additional cut is that we are assuming all along 
LB as a working hypothesis. This predicts, for ideal data, that the velocity vectors reconstructed from the clusters' 
distribution should be aligned with the velocity vectors of mass deduced by POTENT. However, the various random 
errors and systematics in both types of real data analysed here cause deviations from this simple picture. In our 
"standard comparison volume" we restrict the comparison only to points where the velocity vectors of the POTENT 
and cluster fields are broadly aligned with each other. Note that the misalignment constraint may, in principle, affect 



our P estimate and therefore it will be dropped in some of the robustness tests performed in § 4.4 



and S 4.5. A 



"standard comparison volume" is defined trough the set of cuts reported in Table |l[ The two cuts a^^ and turned 
out to be ineffective, in the sense that they are redundant for reasonable choices of the other parameters, and were 
not implemented. The as^ constraint is obtained by scaling as^ by the /3c of BP96. For the sake of consistency, we 
perform all tests with the mock catalogs using the same standard cuts, even though the AO cut does not affect the 
S-5 comparison. The standard volume, Vst, depends on the smoothing applied. For a Gaussian filter of 15 h~^Mpchas 
an effective radius Re ^ 38/i~^Mpc (defined by (47r/3)i?e — Vst)- The rms of ag^ there is ^ 0.18 and of a^^ it is 
~ 150kms~^. Finally, note that the misalignment criterion depends on the the particular velocity field of the generic 
mock catalog and therefore the same standard cuts define slightly different comparison volumes in each of the mock 
Mark III and cluster catalogs tested. 

In what follows we will present the results as the average of the individual results obtained for each catalog, and 
for illustrative purposes also show the results obtained for the mean fields, averaged over the mock catalogs. (The 
volumes corresponding to the individual mock catalogs typically share 80% of their points with the volume defined 
for the average fields). Part of these errors are systematic. The systematic errors can be evaluated by inspecting the 
average of results over the mock catalogs or by comparing directly the average POTENT density and velocity fields, 
to the underlying smoothed fields of the simulation. The top panel of Figure 1 shows this comparison for the G15 
smoothed density fields, at the points of a uniform grid inside the standard volume. The residuals in this scatter plot 
{{Sp) vs. 5t) are the local systematic errors. Their rms value over the standard volume is 0.08. The corresponding 
rms of the random errors {Sp vs. (Sp)) is 0.16. The systematic and random errors add in quadrature to give the total 
error {Sp vs. St), whose rms over the realizations at each point, (jp, is used in the analysis below. 

To quantify the effect of these errors on the determination of /3 we perform a regression of Sp on St, for each 
mock catalog and for the average field, by minimizing the following x^'- 

^2^yi {5p,.-A-BSt,.)' (2) 

The figure shows no considerable systematic deviations from the Sp = St line (the slope of the regression for the 
average field comes out 1.01). The average of the slopes over the 20 mock realizations comes out slightly deviant from 
unity, 1.06, with a standard deviation of 0.17. For the other smoothings G12 and G20, in the standard comparison 



volume, the average of slopes is 1.06 ±0.09 and 0.98 ±0.28, respectively. For other choices of the cuts (see § i.4 § 4.5 
for some illustrative examples) the slope changes by a couple of percent, e.g., for G15, the typical deviation from 
unity is by up to ~ 5% going either way. Finally, only negligible zero point offsets have been detected in all of the 
above comparisons. These results indicate that the final systematic errors are hardly correlated with the signal or 
with themselves, contributing no significant bias in comparisons with density from redshift surveys. 

An analogous comparison has been performed between the supergalactic Y components of the velocity fields at 
the same points. We limit our analysis to this component only as it is the one least affected by the uncertainties of 
the mass distribution near the galactic plane in the forthcoming comparison with cluster velocities. Using the two 
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remaining Cartesian components would require an estimate of possible systematic errors that are uncertain for the 
cluster case (see § 3.3.2). The bottom panel of Figure 1 shows the corresponding scatter diagram of the average 
field vs. the underlying one, again with G15 smoothing in the standard volume. The rms value of the residuals over 
the comparison volume which represent the local systematic error, {{vp) vs. vt), is 75kms^^. The rms value of the 
random errors around the average, {vp vs. { Vp )), is 130kms~^. 

Visual inspection of the figure shows clear signs of systematic errors. The peculiar morphology in the velocity- 
velocity scatterplot reflects correlated velocities within individual cosmic structures. Indeed, the coherence length of 
the velocity field is much larger than for the density field leading to oversampling and correlations among the errors. 
The overall effect on the slope is a bias toward smaller values than unity. The slope of the best-fitting line for the 
average field is in this case 0.81, and the average of slopes over the mock catalogs reflects this as well, giving an 
average slope of 0.80 ± 0.18. The average slope is 0.93 ± 0.13 for the G12 case, and it is 0.76 ± 0.31 for G20. When 
varying the volume, in the G15 case, the bias is typically 12 — 22%. Thus, the velocity comparisons tend, in general, 
to be less robust than the density comparisons and also more sensitive to the smoothing scale. This is probably 
partly due to the larger cosmic scatter in the velocity field, to larger systematic biases in the POTENT analysis (in 
particular the window bias and the sampling-gradient bias) which become more severe for large smoothing scales, 
and to correlation among the errors. 

The POTENT output of the real Mark III data is similarly provided, for the three different smoothings, on a 
Cartesian grid of spacing 5 /i~^Mpc, within a volume of radius 80 h^^Mpc. The errors at each grid point {ag^ and a^^) 
are taken to be the error estimates of the mock catalogs detailed above, i.e., the rms difference over the realizations 
between the recovered fields and the true underlying one. 



3 RECONSTRUCTING THE CLUSTER DENSITY AND VELOCITY FIELDS 

The present analysis is based on the real space cluster distribution and peculiar velocities recovered from the observed 
distribution of Abell/ACO clusters in redshift space. The details of our reconstruction method are described in BP96. 
Here we brieffy describe the data and sketch the main features of the procedure including the error analysis. It is 
worth noticing that in the present comparison with POTENT we are mainly interested in a local region of radius 
~ 70/i~^Mpc, where the reconstruction technique is more reliable than at larger distances. 

3.1 Cluster Data 

The cluster sample used in BP96 contains all the Abell and AGO clusters (Abell 1958; Abell, Corwin and Olowin 
1989) of richness class R > within 250 /i~^Mpc, |6| > 13° and mio < 17 (where mio is the magnitude of the 
tenth brightest cluster galaxy as corrected in PV91). The Abell and AGO catalogs were unified into a statistically 
homogeneous whole-sky sample of clusters using the distance-dependent weighting scheme of P91. The sample used 
in this work contains the same ~ 500 clusters, for which 96% now have measured redshifts, most recently from the 
ESO Nearby Abell Cluster Survey (Katgert et al. 1996, EN AGS). For the remaining ~ 20 clusters the redshifts were 
estimated from the m\o-z relation calibrated as in PV91. The results from this improved sample turn out to be fully 
consistent with the original ones of BP96. 

3.2 Reconstruction of Uniform Cluster Catalogs in Real Space 

Our reconstruction procedure is of two steps. First, Monte Carlo techniques are used to correct for observational 
biases and return a whole-sky distribution of clusters in redshift space. Then, this distribution is fed into an iterative 
reconstruction procedure (similar in spirit to Yahil et al. 1991) which assumes linear GI-I-LB to recover the real-space 
positions and peculiar velocities of the clusters. 

The main observational biases arise from a systematic mismatch between the Abell and AGO catalogs and from 
the latitude-dependent Galactic obscuration; the radial selection is not an issue because it is quite uniform in the 
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volume relevant for our analysis. To minimize the possible systematic errors in the model cluster velocity field we 
need to unify the Abell and AGO catalogs into a statistically homogeneous whole-sky sample of clusters. We obtain 
this by using the distance-dependent weighting scheme of PV9f which enforces the same number density in equal 
volume shells for the two cluster populations. The number of radial shells is left as a free parameter. To correct 
for Galactic obscuration, we generate a set of cluster catalogs of uniform sky coverage, by adding a population of 
synthetic clusters. The Galactic obscuration at |6| > 20° is modelled by a cosecant law, 'P{h). As in BP96 we have 
chosen two different sets of absorption coefficients to account for observational uncertainties. Synthetic clusters are 
added with a probability proportional to V{b) such that they are spatially correlated with the real clusters according 
to the observed cluster-cluster correlation function. Within the ZoA, |6| < 20°, the volume is filled with synthetic 
clusters, in bins of redshift and longitude, by cloning the cluster distribution in the adjacent latitude strips outside 
the ZoA. Both real and synthetic, Monte Carlo generated clusters are mass weighted to determine the density field 
to be used in equation (^) below. The mass of each real cluster is proportional to the number of galaxies per cluster 
listed in the Abell catalog. The mass of synthetic clusters is set equal to the mass of the real ones. 

The redshift space distortions are corrected by an iterative procedure based on linear theory and linear biasing. 
Equation (|l|) can be inverted to yield 



This is used, in each iteration, to compute the radial peculiar velocities of the individual clusters, u, in the LG frame, 
and correct for their real distances, r, via r = cz — u. To avoid strong nonlinear efi'ects, the force field generated 
by the point-mass clusters is smoothed by a top-hat window of radius 15/i^^Mpc, chosen to be comparable to the 
cluster-cluster correlation length (see BP96). A meaningful comparison with POTENT, in view of the large mean 
separation between clusters (~ 25/i~^Mpc), requires that we smooth further the density and velocity fields. As input 
to this procedure one has to assume a value for Pc, which affects the peculiar velocities but has only a weak effect 
on the real-space distance as long as /3c is in the right ballpark (see BP96). We have assumed — 0.21 based on 
matching the dipoles of the GMB and the cluster distribution (e.;;., BP96). 

3.3 Errors in the Cluster Positions and Velocities 

Ideally, we would like to implement the same POTENT error assignment procedure also for the cluster case. However, 
intrinsic difficulties in modelling the Abell/ACO selections criteria hamper the compilation of mock cluster catalogs 
from the K96 simulation. A problem which is made worse by the size of the N-body computational volume which is 
smaller than the one spanned by the real clusters' population. We choose instead to evaluate cluster errors using a 
hybrid scheme in which: 

• The Monte Carlo procedure of adding synthetic clusters and varying the parameters is also used to estimate the 
random and systematic errors that arise both from the uncertainties in modelling the observational biases and the 
approximations in the reconstruction. 

• A mock catalog analysis, similar to the one used to assess the POTENT errors, is implemented to quantify the 
additional random errors that arise from the sparseness of the clusters' sampling. 

3.3.1 Monte Carlo Analysis 

The reconstruction procedure depends on a number of parameters that are only weakly constrained by observational 
data or theoretical arguments, such as the galactic absorption coefficients, the force smoothing length, and the 
weighting scheme used to homogenise the Abell and AGO catalogs. BP96 evaluated the sensitivity of the derived 
density and velocity fields to these parameters by allowing the parameters to vary about the standard set defined in 



The total uncertainties in the cluster positions and in the radial velocities, estimated in the CMB frame, arise 
from several different sources: 




(3) 



their Table 2. 
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• Intrinsic errors of the reconstruction procedure, which we estimate by the standard deviation of the cluster 
distances over 10 Monte-Carlo realizations of the same choice of parameters. 

• Observational errors, accounting for the freedom in the values of the free parameters, which we estimate by the 
standard deviation of the distances over reconstructions with different sets of values for the parameters. 

• Shot-noise error (eq. [20] of BP96), due to the uncertainty in the mass per cluster which we assume proportional 
to the number of galaxies listed in the Abell catalog (see BP96). This error is estimated to be of ~ 70kms~^ (one 



dimensional). The shot-noise due to the sparseness of the mass tracers will be estimated numerically in § B.3.2. 

• Weight uncertainty. This error accounts for uncertainties in the relative weighting of Abell versus AGO clusters 
to correct for systematic differences between these catalogs. For this purpose, we have performed 10 different recon- 
structions in which the weights were randomly scattered about the standard weights of BP96 (their eq. [5]), following 
a Gaussian distribution of width that equals the Poisson error in the relative number densities of Abell and AGO 
clusters at every given distance. The estimated typical weight uncertainty turns out to be ~ 85kms^^. 

• Projection uncertainty. A worry when using the Abell/AGO clusters for statistical purposes is the contamination 
of cluster richness due to projection of foreground and background galaxies (e.g., Dekel et al. 1989 and references 
therein). The resulting uncertainty in the galaxy count per cluster has been recently estimated (Van Haarlem et al. 
1997, using N-body simulations; Mazure et al. 1996 using the ENAGS survey) to be ^ 17%. To translate this error 
into a distance error, we have performed 10 reconstructions where the cluster richness were randomly perturbed by 
a 17% Gaussian, yielding an error of ^ SOkms"'^. 

An upper bound to the total error for each Abell/ AGO cluster can be estimated by adding in quadrature all 
the above errors, as if they are all independent. This results in an average error of 254 km s^^, with a large spread 
of ±117kms~^. If add in quadrature only the observational, intrinsic and shot-noise errors, which are independent 
of each other, the average error drops only to 218 km s"'^, indicating that our upper bound is not a far over-estimate 
of the true error. Figure 2 shows the distribution of the total reconstruction errors in the line of sight component of 
the peculiar velocities for the ~ 500 clusters of our subsample. Unlike the intrinsic ones, observational, shot-noise, 
weighting and projection errors are isotropic. Therefore, they are also representative of the uncertainties along the 
supergalactic Y component of the velocity fields and they will be used to estimate the cluster velocity errors (t„c in 



3.4) 



For the comparison with POTENT we should identify the regions where the reconstruction from clusters is 
reliable. The errors are naturally larger in regions where the fraction of observed clusters is lower. This effect is 
clearly seen when we plot the intrinsic error per cluster as a function of Galactic latitude (Fig. 3). As expected, no 
radial dependence has been detected for the errors within the volume used for the present analysis. 



3.3.2 Mock Catalog Analysis 

Due to their low number density, clusters of galaxies sparsely sample the underlying density and velocity fields. This 
introduces an intrinsic scatter, sometimes termed 'shot noise' as well, when comparing the cluster and the mass fields. 
This is closely related to the expected scatter from the stochasticity in the bias relation (e.g. Dekel & Lahav 1999). 
This sort of random errors are not included a-priori in our cluster error estimates, but need to be accounted for in 
the actual comparisons of the clusters fields to the POTENT reconstructions. 

We further assess the reliability of the cluster fields, and get a crude estimate of the additional scatter, using 
the same N-body simulation that was the basis for the Mark III mock catalogs. In the present case, however, we 
obtain just one mock catalog of clusters from the simulation. We use a friends-of-friends algorithm to identify groups 
among the particles in the simulation. The richest groups above some threshold, fixed so as to have the same number 
density as the Abell/AGO clusters, are identified as the mock clusters. To mimic the properties of the real cluster 
distribution we need to extend our mock sample out to 250/i^^Mpc. Since this exceeds the size of the simulation, 
we obtain the mock cluster distribution by duplicating the clusters within the computational box using the periodic 
boundary conditions. The peculiar velocities from clusters are then computed from equation (^) and both the cluster 
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density and the velocity fields are smoothed at the points of a cubic grid with 5 /i^^Mpc spacing. The smoothed cluster 
fields are then compared with the true underlying fields of the simulation, smoothed on the same scale. Under the 
GI+LB assumptions, the fields are simply related by the biasing factor between clusters and mass in the simulation, 
both in the density case and for the velocities (for our = 1 simulation). We use the same "standard comparison 
volume" considered for the POTENT vs. true comparisons and used later for the POTENT vs. cluster comparisons. 
Figure 4 shows the results for the G15 case. 

The slopes of the best fitting lines for the 5-5 comparison (0.32, top panel) and for the Vy-Vy comparison (0.30, 
bottom panel) have been estimated by assigning equal weight to all points in the plots. They are a measure of the 
relative biasing between clusters and mass {hc~^ , when regressing clusters on true fields). 

The average value of the 20 volumes defined by the POTENT mock catalogs with the same criteria is 0.31 ± 0.03 
for the density fields and 0.33 ± 0.05 for the velocities. For general variations in the comparison volume, values of 
0.30 — 0.38 are typically obtained, with a slight tendency of the values obtained from velocities to be higher than those 



obtained from densities, within this range. Unlike the POTENT analysis in § 2.1, we do not know the 'true' expected 
slope for the mock clusters and this sort of comparisons are in practice a way to define it. Variations between the 
values obtained from densities and from velocities may arise because of the larger cosmic scatter for the velocities and 
uncertainties in modelling the cluster distribution outside the computational volume. Another possible cause for the 
mismatch is the already mentioned strong correlation among the errors in the Vy-Vy analysis. Note that the difference 
between the slopes obtained from the 5-5 and the Vy-Vy is smaller than the outcome of the POTENT analysis in § 2.1, 
meaning the the various error sources affecting the Vy-Vy comparison tend to compensate each other. 

The distribution of the distant clusters may significantly affect the cluster velocities while it is almost irrelevant 
when computing the smoothed density field within 70 Mpc. We therefore regard hc~^ ~ 0.31 as the 'true' value 
for the G15 standard case. Similar values are obtained for hc~^ with the other smoothing scales. 

A considerable scatter about the regression lines is found in both the 5-5 and the Vy-Vy comparisons. Since 
the clusters in this case are free from the observational and modelling errors, this scatter is a manifestation of the 
additional inherent scatter in the cluster fields mentioned above. For G15, the detected scatter is ct^'"' = 0.36 in the 
density case and (Ti,*"* — 300kms~^ for the velocities. The change with smoothing scale is as expected: for G12 the 
scatter is larger (0.46 and 400, for densities and velocities, respectively) and for G20 the scatter is smaller (0.24 and 
200). These estimates are quite robust to changes in the comparison volume. 

In what follows, we adopt these dispersions as a measure of the intrinsic scatter of the cluster fields, for the 
POTENT-cluster mock tests in § 4.2, and also for the comparisons with the real data. The drawback of the latter 
assumption is the fact that the mock clusters do not accurately match the Abell/ACO cluster distribution. Most of 
the mock clusters do not correspond on a one-to-one basis to the Abell/ACO ones, and are less spatially correlated. 
Furthermore, since the mock clusters are identified from a simulation based on IRAS galaxies, which tend to avoid 
high density regions, they might represent density peaks of a systematically lower amplitude with respect to those of 
the Abell/ACO clusters. Finally, the a-similarity could be even more severe for the velocity calculations because of the 
duplication procedure adopted outside the N-body computational volume. Still, not having a better way to accurately 
determine this scatter for the real data, and understanding that its magnitude mainly depends on the sparseness of 
clusters and smoothing adopted, we believe that our approach does give a crude estimate of the effect. It is interesting 
to check the plausibility our results by comparing them to the analytic estimates of shot noise computed according 
to Yahil et al. (1991). For G15, and within a radius of 70 h^^Mpc, we obtain (Ja""' = 0.42 for the 5-5 case and 
o"t,'"* = 1540/3c kms~^ for thevy-Vy one. Scaling to (5c = 0.21 value of BP96 we obtain that in both cases the analytic 
shot noise is close, although somewhat larger, to the scatter in the simulations. The explicit assumption made here is 
that the intrinsic scatter found in the mock simulation is representative for the real clusters too, and is independent 
of the other sources of error and the underlying field. The plausibility of these hypotheses will be assessed a posteriori 
when comparing the real cluster and POTENT fields. 



As for the POTENT analysis (§ 2.1), the cluster velocity analysis has been limited to the supergalactic Y 



component, that is less prune to systematics. More extended mock catalog analyses, based however on the distribution 
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of IRAS galaxies, have demonstrated this point (Branchini et al. 1999). The other two Cartesian components are 
affected by systematic errors, arising from the cloning procedure that is used to fill the ZoA. Given the larger extent 
of the ZoA in the cluster case, we expect comparable, if not larger, systematics to affect the cluster velocity field. 
Correcting for this bias would require an error analysis based on more realistic Abell/ACO mock catalogs that are 
not currently available. Therefore, we restrict our analysis to a Vy-Vy comparison, under the working hypothesis that 
the Y component of the cluster velocity field is only affected by random errors. 



3.4 Smoothed Density and Velocity Fields 

For the purpose of the comparison with POTENT, we compute smoothed density and velocity fields at the points of 
a cubic grid with spacing 5 /i~^Mpc inside a box of side 320 h~^Mpc centered on the Local Group. 

We first generate 20 Monte-Carlo realizations of our standard model as described above. To mimic the effect of 
observational errors and shot noise, we perturbed the cluster distances with a Gaussian noise of 150 kms-^ This 
value is slightly larger than the sum in quadrature of average observational and shot noise errors and corrects for 
the positive tail in the error distribution similar to the one observed in Figure 2. The intrinsic error is not included 
here because it will enter implicitly when we later average over the 20 fields. A Cloud-in-Cell (CIC) scheme is used 
to translate the discrete cluster distribution into a density field at the grid points. The peculiar velocities at the grid 
points are recomputed from the reconstructed cluster distribution via linear GI+LB, with the force field smoothed 
by a small-scale top-hat window of radius 5 /i~^Mpc. We minimize the scale of the top-hat force smoothing in order 
to eventually end up as close as possible to Gaussian smoothing on scales > 12 h~^Mpc. 

The 20 fields are then smoothed further with a Gaussian window of a larger radius {Rg ~ Ri)'^'^ , where Ri is the 
Gaussian smoothing radius equivalent in volume to the 5/i~^Mpc top-hat force smoothing. As mentioned already, 
we try three different smoothing scales, of Rs = 12, 15 and 20/i~^Mpc. These 20 fields are averaged to give the final 
smoothed fields used in the next section. Each of these 20 fields is affected by observational and shot noise errors. 
Moreover, since they represent 20 different Monte-Carlo realizations of the cluster distribution for the same choice 
of parameters, we can account for the intrinsic errors by the very same averaging procedure. Indeed, the standard 
deviation over the 20 realizations represent the cumulative effect of intrinsic, observational and shot noise errors 
The only contribution left to the total error budget is the intrinsic scatter estimated in § j.3.2| . 



discussed in 



3.3.1 



This is modeled as a Gaussian noise and is added in quadrature at each gridpoint. Therefore, the error estimates 
for the smoothed cluster fields, as^ and a^^ , are obtained by taking the standard deviation over the 20 catalogs and 
adding in quadrature a Gaussian noise of amplitude equal to the intrinsic scatter. 



4 MEASURING f3c 

We determine /3c in two different ways, via 5-5 and Vy-Vy comparisons. Because the 5-5 comparison is local, it avoids 
the incomplete sky coverage in the ZoA, but it uses only a small number of clusters (18 within 70/i^^Mpc). The 
Vy-Vy comparison, on the other hand, involves an integral over the cluster distribution in an extended volume, but it 
suffers from a large uncertainty due to the unknown cluster distribution in the ZoA and beyond the sample's edge. 
Therefore, the two methods are expected to suffer from different biases and can provide us with two estimates of (3^ 
that are somewhat complementary. 

We wish to restrict the quantitative comparison to the regions in space where both the errors in the cluster and 
POTENT fields are reasonably small. On the other hand, we wish to maximize the number of independent volumes 
compared in order to minimize the cosmic scatter. We therefore need to optimize our choice of comparison volume, 
and test the robustness of the results to changes in this volume. The comparison volume has already been introduced 
The natural parameters for defining it the are the measure of poor sampling of the Mark III catalog, Ra, our 



2.1 



estimate of the random errors in the POTENT and cluster density fields {os^ and as^ ) and the corresponding errors 
in the velocity fields {(Jvp and (Ti,^). The latter turned out to be ineffective cuts, in the sense that they are redundant 
for reasonable choices of the other parameters. Our 'standard' cuts according to these parameters are reported in 
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Table As mentioned already in § 2.1 we impose as well a constraint on the misalignment angle between the cluster 
and POTENT velocity vectors. This serves as an additional classifier of "good" points for the comparison, and helps 
avoiding regions where we might have large, perhaps unaccounted for, errors. In our main analysis we restrict the 
comparison to points with a maximal misalignment angle of AS < 45° . We later relax this constraint and verify the 
robustness of the results. We take the G15 smoothing as our standard case, with the above set of criteria defining 
our 'standard' comparison volume. 

4.1 The Pc Fitting Method 

The assumption underlying the /3c estimations is that the density and velocity fields recovered above are consistent 
with the model of GI+LB. In this framework the POTENT and the cluster fields are linearly related: 

p = /3cc + A, (4) 

where p and c stand for the POTENT and cluster and represent either 5 or the supergalactic Y velocity component. 
The cluster errors, Oc, are comparable to the POTENT errors, (jp. The best-fit parameters are therefore obtained by 
minimising the quantity 



2 ^ {P^-A- l3,Ci) 

i— 1 

where the subscript i refers to any of the Ntot gridpoints within the comparison volume. Since the fields have been 
smoothed on scales much larger then the grid separation, these points are, however, not independent. As in Hudson 
et al. (1995; see also Dekel et al. 1993), we estimate the effective number of independent points, Neff, as 

Ntat Ntat 

N^/f = N^ot J2 J2 ^M~rl/2Rl) , (6) 

where Vij is the separation between gridpoints i and j. This expression weighs the dependent grid points taking 
into account properly the finite comparison volume and its specific shape. This estimate is thus more accurate than 
the simplistic ratio of the comparison volume over the effective volume of the smoothing window, which assumes an 
infinite comparison volume. We account for the oversampling problem by using an effective statistics defined by 
Xe// = {Neff /Ntot)x^ J which is equivalent to multiplying the individual errors by the square root of the over-sampling 
ratio Ntot /Neff. The assumption we make is that this new statistics is approximately distributed like a with N^ff 
degrees of freedom. In what follows we use it to assess the errors in f3c and A. 

4.2 Testing the Comparison 

Before performing the comparisons with the real data, we wish to further quantify the possible systematics that might 
enter, verify the validity of the smoothing scheme adopted and find the optimal smoothing scale for the comparisons. 
Also, we would like to understand whether the intrinsic differences between the density and velocity field comparisons 



can affect the results. We do this by comparing the mock POTENT fields of § 2.1 with the mock cluster ones from 



§ 3.3.2. The results of the density and velocity comparisons, within the standard comparison volume, for each of our 
three smoothing scales, are reported in Table ^. The values quoted in the table are all the mean values averaged over 
the results of the 20 mock POTENT catalogs. The G15 case is illustrated in Figure 5, where the average POTENT 
fields are compared to the cluster ones. 



The values obtained for /3c in the G15 case are encouragingly close to the "true" value of 0.31 (§ 3.3.2). The 
typical range of values obtained when altering the chosen volume for the comparison, from densities and velocities, 
lie in the range 0.26 — 0.35. 

There are small differences for the other smoothings, with a tendency towards smaller /3c for the G20 case. All 
variations are, however, well within the formal Xeff error-bars. In all cases, the zero-point found is hardly significant 
and is consistent with zero, within the error-bars. 
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The /Neff values (defined as S in tlie table) are only slightly smaller than unity for the density comparisons, 
well within the accepted range given the small number of degrees of freedom (1±0.4 for G15), but they are significantly 
smaller for the velocities. A similar trend was also found when comparing the POTENT fields with the true N-body 
ones. This may indicate that model velocities are more correlated than the scaling to A^e// obtained for the densities. 
If true, then our Xe// distribution deviates from that of a with Neff degrees of freedom, resulting in an error 
overestimate and in a /Neff < 1 for velocities. Note that this effect is also associated with the application of the 
alignment constraint, and without it the value of x^ increases somewhat. So, perhaps for aligned vectors our errors 
are over estimated. Another, possibly not exclusive, explanation, which is also suggested by the POTENT vs. N-body 
fields comparison, is the existence of systematic errors that do not average to zero. 

The formal error-bars obtained for /3c, in all cases, are significantly larger (by a factor of ~ 2) than the spread of 
values obtained from the 20 Mark III mock catalogs. Again, this may suggest that our effective x^ statistic recover 
the correct slope but overestimates the errors on /?c also for the 5-5 comparison. 

In summary, the important conclusion from the comparison of the mock data is that our method provides a fairly 
reliable estimate of /3c with no gross biases. The results are fairly robust to changes in the comparison volume and 
the smoothing scale. The reasonable values obtained for the densities versus the too-low values for the velocities, 
incline us to regard, also here, the density comparison as the more reliable one. 

Our purpose is to constrain /3c in a meaningful way by comparing the density and velocity fields extracted from 
the cluster distribution with the fields recovered by POTENT from peculiar velocities. The competing obstacles are 
the very sparse sampling of the underlying density field by the clusters, on the one hand, and the limited volume 
sampled by peculiar velocities, on the other. The former dictates the use of a large smoothing scale, because small- 
scale structure is not traced properly by the clusters, while the latter calls for a relatively small smoothing scale, in 
order to minimize the cosmic scatter associated with the number of independent volumes, and the systematic biases 
in the method. The rough agreement between the results of diflferent smoothings is encouraging. Large smoothing 
scales, such as the G20 case, are perhaps more pruned to systematics, and in any case, since they reduce the number 
of independent data points, they can only constrain /3c very weakly with large uncertainties. These considerations, 
along with the fact that it matches the intercluster separation, made us choose the G15 filter as our standard for the 
real data analysis. As mentioned before, the formal errorbars from the Xe// statistics may overestimate the actual 
uncertainty in the results, but we conservatively choose to stick with these. 

4.3 Visual Comparison of Maps 

Figures 6 display the G15 density and velocity fields (in the CMB frame) from the clusters (left) and POTENT 
(right) reconstructions of the real data in three slices parallel to the Supergalactic plane, within a sphere of radius 
80/i~^Mpc about the Local Group. The clusters' densities and velocities are scaled by j3c = 0.21. The heavy line 
delineates our standard comparison volume. 

The similarity between the two density fields is evident in most regions. In both fields, the dominant features axe 
the Great Attractor (on the left), the Perseus-Pisces supercluster (on the right), and the great void in between. On 
the other hand, the Coma supercluster, seen in the clusters map near {X, Y) « (0, 70), is not reproduced at the same 
position in the POTENT map. Differences are also seen in the upper-right quadrant of the Z = — 25ft.~^Mpc plane. 

There is also some qualitative agreement between the velocity fields, but it is less striking. The main features 
common to the two fields are the convergences into the Great Attractor and into Perseus Pisces. The main difference 
is an additional bulk flow from right to left for the POTENT field, apparent in the three slices. Another feature absent 
in the POTENT field is the infall into Coma seen in the cluster field. Note that the main discrepancies lie outside 
the comparison volume, in regions where the errors are expected to be large at least in one of the reconstructions. 
These regions will be excluded from the quantitative comparison below. 

It is also worth noticing that the density-velocity maps for the clusters are very similar to those obtained by 
Scaramella (1995b) from the same Abell/ACO cluster catalogues but using a somewhat different technique. 
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4.4 Estimating fic by a Density Comparison 

We perform the 5 — 5 comparison within the standard volume. The errors in the POTENT field have been evaluated 



2.1 



and for clusters we use the error estimates of § 3.3. The results for the three smoothing radii are displayed 
in the first three rows of Table For the preferred G15 case we find /3c = 0.20 ± 0.07, and the best-fit value stays 
essentially the same for the other cases (with the errorbar increasing with the smoothing, due to the smaller number 
of effective independent points). No significant zero point offset is found in any of the cases. The 5-5 scatterplot is 
displayed in Figure 7 for the G15 case. The solid line is the best-fit from the minimization. 

We have tried several variants of the comparison volume, in order to check the sensitivity of our results. Two 
representative examples are reported in the last two rows of Table ^ In the forth case there we have considered the 
original standard volume but with a stricter Ri cut, R4, < 10^~^Mpc. The last column shows the results for the 
most interesting experiment, i.e. the one in which the misalignment constraint has been removed. The results of these 
tests all confirm the robustness of the /3c — 0.2 value. For the density comparisons, we generally get xtf f / ^ef f ~ 1, 
indicating a good fit. Note, that relaxing the constraint on the misalignment angle more than doubles the number of 
gridpoints considered. 
As outlined in 



3.3.2, the present results have been obtained assuming that the scatter found in the mock cluster 



fields is representative of the intrinsic scatter in the real case and that it is independent of the other sources of error 
which form as^. The resulting Xe// values are an indication that these are indeed fair assumptions. 



4.5 Estimating /3c by a Velocity Comparison 

As already pointed out, it is important to perform the POTENT-cluster regression for the velocities on the ground 



of its complementarity with the 5 analysis. Also, as we have already discussed in § 3.3.2 , we limit the comparison to 
the supergalactic Y component which is the more robust of the components. We use the same minimising procedure 
adopted in for the 5-5 comparison. The results are displayed in the right half of Table ^ For our standard G15 case 
the result is now /3c = 0.25 ± 0.05, somewhat higher than the density case, but still consistent within the errorbars. 
The scatterplot for this case is shown in Figure 8. 

As was the case for the 5-5 comparisons, no significant offset is detected, and the /3c value is quite robust for the 
different smoothing scales and under variations of the comparison volume (the changes in the resulting /3 are well 
below the la significance level). Note again the peculiar morphology in the scatterplot, arising from the coherency 
in the peculiar velocities within independent cosmic structure. Although it may seem that the POTENT and cluster 
velocity fields differ in a large scale bulk flow component, more quantitative, volume limited comparisons performed 
by Branchini Plionis and Sciama (1996) and Branchini et al. (1999) have shown that the two bulk flows agree in 
amplitude and direction for a value of /3c ~ 0.21 

It is especially interesting to check the effect of removing the alignment constraint (the fifth case in the table). 
This is a demanding robustness check since it extends the Vy-Vy comparison volume to points for which the velocity 
vectors can be severely misaligned. It is encouraging that, even in this case, the slope of the best fitting line changes 
only by 4 %. The /N^ff values lie somewhat below unity for all the cases explored but for the one in which we have 
removed the alignment constraint. In this las case we obtain /Neff ~ 1 both for the 5-5 and Vy-Vy comparisons. A 



similar behavior was also obtained for the mock comparisons (§ 4.2). 

The errors cr^^ obtained from the real analysis are smaller than those obtained from the mock and listed in 
Table ^ Even accounting for the difference in the values of /3c, the two error estimates differ by a factor of ~ 2. This 
mismatch probably arises from the characteristics of the mock velocity fields. Indeed, the small computational box 
used in the original K96 simulation and the constraint of having a vanishing bulk velocity on the scale of the box, 
produce a remarkably quiet velocity field with a bulk velocity of only 100 kms^^already on a scale of 40 h~^Mpc. 
This velocity field has been used to estimate the POTENT and part of the cluster velocity errors. Real velocities, 
however, are larger than the mock ones and these uncertainties probably underestimates the errors for the real case, 
leading to the smaller ct^^ value listed in Table ^ 
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The Vy-Vy comparison described above has been performed in the CMB reference frame. Predicting velocities 
from galaxy redshift surveys is commonly done though in the LG frame, in order to minimize the influence of mass 
concentrations from outside the sample volume. The LG frame might therefore be considered the natural frame 
in which to perform comparisons with reconstructed velocities. In our case, the velocities arc reconstructed from 
the far-extending cluster catalog, which alleviates the above problem and wc regard a CMB comparison as reliable. 
Furthermore, performing the comparison in the LG frame would introduce extra complexities, requiring a somewhat 
ad-hoc transformation to a common LG frame for both the cluster and POTENT velocity fields. As a crude test 
of the sensitivity of our results to changes in the framework of reference, we shift both velocity fields to the cluster 
LG frame, as defined by the smoothed cluster velocity at the origin (with a reasonable choice for /3c)- Alternatively 
we consider the peculiar velocities relative to the central observer of each reconstructed velocity field independently. 
The standard G15 comparison of the Y components gives /3c = 0.24 ± 0.05 and 0.25 ± 0.04 respectively, for these two 
cases, demonstrating once more the robustness of our result. 



5 CONCLUSIONS 

We have used the smooth matter fiuctuation field obtained by applying the POTENT machinery to the Mark 111 

dataset and compared it to the density field deduced from the Abell/ACO cluster distribution. A similar comparison 
has also been performed between the reconstructed cluster velocities and those from the Mark 111 catalog, smoothed 
on the same scale. We have performed a careful error analysis using mock galaxy and cluster catalogs derived from 
N-body simulations. The mock catalogs used in our POTENT error analysis were especially designed to reproduce the 
Mark 111 characteristics. Uncertainties in the cluster fields, on the other hand, were evaluated using a hybrid procedure 
which extends the Monte-Carlo error analysis of BP96 and is complemented with a similar to the POTENT mock 
catalog analysis. 

Cluster and POTENT fields show remarkable similarities within 70/i^^Mpc, while their major discrepancies 
are usually confined into regions where the cluster or the POTENT reconstructions are known to be unreliable. 
Quantitative comparisons between cluster and POTENT fields have been performed in an attempt to estimate the 
cluster /9 parameter. The results arc quite robust and for the standard G15 case wc find Pc ~ 0.20 ± 0.07 from the 
S-S regression, and a somewhat larger value of f3c = 0.25 ± 0.05 from the Vy-Vy case. This systematic discrepancy 
is within the la significance level, but it is present in all the cases explored. We therefore choose to quote a joint 
estimate for /3c of 0.22 ± 0.08. 

Some differences between the two values are not unexpected given the different nature of the comparisons. A 
similar regression based on the mock catalogs showed that some discrepancies do exist. However, in the mock tests 
the difference between the two values was of smaller magnitude and in the opposite direction. The different trends 
between the real and mock results could arise from the different modelling of the mass distribution outside the 
sampled regions, which can affect the cluster velocity field. There are other indications for regarding the 6 — 6 results 
as being more reliable. The x^/^eff values for the density comparison were around unity, while systematically lower 
values were obtained for the velocities. Also the POTENT velocity field was found in the mock catalog analysis to 
suffer of more biases. 

The present analysis suggests a value of /3c — 0.20 — 0.25 for clusters, in accordance with previous estimates. 
The distribution of clusters is expected to be biased with respect to the distribution of galaxies with a biasing factor 
6cg — 3 — 4 {e.g., from the different correlation lengths obtained for clusters and for galaxies; Bahcall & Soneira 1983, 
Huchra et al. 1990). Peacock and Dodds (1994) find such values for the biasing factors, derived from the ratios of power 
spectra calculated for different datasets. Their quoted relative biasing factors for Abell clusters, radio galaxies, optical 
galaxies and IRAS galaxies is 4.5 : 1.9 : 1.3 : 1, respectively. Recent results from a comparison of the cluster density 
and velocity fields to the fields recovered from the PSC« redshift survey constrain this parameter to beg = 4.4 ± 0.6, 
with respect to IRAS galaxies (Branchini et al. 1999). Joined with our constraint on /3c this implies /3/ ~ 1 with, 
however, a 1-cr uncertainty of ~ 50 %. Although our analysis cannot provide us with a firm /3c determination, due to 
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the large uncertainties associated with the 5-& and Vy-Vy comparisons, it leads toward a value of /3c which is consistent 
with an Einstein de Sitter universe for a reasonable cluster linear bias parameter of he ~ 4.5. Our value of be is a 
linear fit to the 5-5 and Vy-Vy scatterplots. Under the assumption of Linear Biasing he represents the relative biasing of 
Abell/ACO clusters with respect to the underlying mass density field. Linear Biasing, however, needs not to be a good 
approximation for clusters since the large value of be causes the LB hypothesis to break down in those regions where 
5 < —be^- Methods to measure the degree of nonlinearity in the biasing relation have been recently developed and 
applied to galaxy distribution (e.g. Lemson et al. 1999 and Narayanan et al. 1999, Sigad, Dekel and Branchini 1999) 
but not yet to clusters of galaxies. However, there are indirect evidences of the small deviation from linear biasing. 
The first one comes from the visual inspection of the 5-5 scatterplot in Figure 7 which does not deviate appreciably 
from LB expectations (the linear fit). A more convincing evidence of the small deviations from LB approximation is 
obtained by performing the 5-5 comparison for various smoothing filters. Increasing the smoothing length decreases 
the amplitude of density fiuctuations and reduces the size of those regions in which the constraint 5 < —b'^ cause 
the LB model to fail. As shown in Table ^, increasing the smoothing radius from 12 to 20/i~^Mpc doesn't change 
Pe significantly showing that regions where LB does not apply play a little role in our analysis. As a consequence, 
and for all practical purposes. Linear Biasing is a good approximation on the scales relevant for our analysis and and 
therefore we can regard be as the biasing parameter for Abell/ACO clusters. 
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Cluster versus POTENT fields 17 



Table 1. The Standard Cuts. Column 1: ag^, Max. error in POTENT S; Column 2: ag^, Max. error in Cluster S; Column 
3: i?4, Max. distance from the 4-th neighbouring object in h^'^Mpc; Column 4: R, Max. radius in h~^Mpc; Column 5: |b|, 
Min. galactic latitude in deg. ; Column 6: A6, Max. misalignment in deg. ; 







i?4 


\b\ 


R Ae 


0.3 


1.43 


13 


20° 


70 45 ° 



Table 2. fSc from the POTENT vs. cluster real data comparisons. The symbols and units are the same as in Table ^. The first 
three rows refer to the standard comparison volume. The last two rows are for G15, when varying the comparison volume (see 
discussion in the text). 



Rs Ntot N^ff Re 13^ a^^ a^^ /3J a^^ A" a\ 5" 

12 1239 11.7 33 0.21 0.06 -0.02 0.06 1.28 0.26 0.05 28 59 0.75 

15 1855 10.6 38 0.20 0.07 -0.04 0.07 1.04 0.25 0.05 27 60 0.69 

20 1537 5.0 36 0.22 0.12 -0.04 0.06 1.23 0.23 0.07 -4 83 0.58 

15 1429 8.5 35 0.19 0.07 -0.04 0.07 1.00 0.26 0.06 22 67 0.73 

15 4286 19.1 50 0.20 0.06 -0.02 0.04 1.01 0.26 0.05 9 48 1.06 



Table 3. fic from the POTENT vs. cluster mocks comparisons. Column 1: the smoothing radius in /i~ Mpc; Column 2: 
Ntot, the number of gridpoints within the volume; Column 3: N^ff, the effective number of independent volumes; Column 
4: the effective radius in h~^Mpc; Column 5: /3^, /3c from Column 6: o"^5, the dispersion; Column 7: , 

zero point offset for 5-5 ; Column 8: ct^, zero point dispersion; Column 9: 5"^ = X^ff/^eff from 5-5; Column 10: /S", /3c 
from Vy-Vy; Column 11: cr^" , the 13^ dispersion; Column 12: A", zero point offset for Vy-Vy\ Column 13: ct^, zero point 
dispersion; Column 14: = Xe///-'^e// from Vy-Vy. 
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FIGURE CAPTIONS 

Figure 1 Systematic errors in the POTENT analysis. The POTENT fields recovered from the noisy and sparsely 
sampled mock data are compared with the "true" G15 fields of the simulation. The comparison is at uniform grid points 
within our "standard comparison volume" of effective radius 40 ft~^Mpc. Plotted, in both cases, is the POTENT field 
averaged over the 20 realizations. Top: The POTENT density field vs. the true density field. Bottom: The POTENT 
supergalactic Y-component of the velocity field vs. the true velocities in the simulation. 

Figure 2 Histograms of the global random+systematic errors from the Monte Carlo analysis. The plot shows 
the frequency of the uncertainties on the line of sight component of the reconstructed cluster velocities. Units are 
kms~^. 

Figure 3 Environmental dependency of Monte Carlo errors. Dependence of the intrinsic errors on galactic 
latitude. Abell/ACO clusters in the northern galactic hemisphere are shown as filled dots while open dots represent 
southern clusters. 

Figure 4 Random errors in the cluster mock analysis. Cluster mock density and velocity fields are obtained from 
the clusters identified in the K96 N-body simulation and compared to the underlying fields, all with G15 smoothing. 
The comparison is at the same points shown in Fig. 1 Top: 5-5 comparison. Bottom: Vy-Vy comparison. The solid 
lines corresponds to the average best-fit lines. The scatter about the fitting lines is an estimate of the intrinsic scatter 
in the cluster fields. 

Figure 5 Comparing the mock POTENT and cluster fields. The averaged mock POTENT fields of 



2.1 



compared with mock cluster ones of § 3.3.2, within the standard comparison volume. Top: 5-5 comparison. Bottom: 
Vy-Vy comparison. The values quoted are the average over the 20 catalogs of /3c and of its estimated error from the 
fit. The solid lines represent the best-fits averages. 

Figure 6 Density fluctuations and projected velocity field in supergalactic X-Y planes. The Mark III-POTENT 
case is shown on the right and the cluster fields on the left, all G15 smoothed. The density contour spacing is 
A(5 = 0.15, solid contours refer to overdense regions while dashed contours refer to negative overdensities. The thick 
line indicates the 5 = contour. The heavy line defines the standard comparison volume. The length of the velocity 
vectors have been drawn on the scale of the plot. The cluster density fluctuations and velocities are scaled by /3c = 0.21. 
Top panel shows the plane deflned by supergalactic Z = -1-2500 kms~^, Middle panel shows the supergalactic plane 
of Z = Okms~^, and the lower panel the plane defined hy Z = — 2500kms~^. 

Figure 7 POTENT versus cluster G15 density field from the real data, at gridpoints within the comparison 
volume. The solid line results from the linear best fit. 

Figure 8 POTENT versus cluster G15 velocity field from the real data. Only the Y supergalactic components 
at gridpoints within the comparison volume are considered. The best-flt line is marked as well. 
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